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Abstract. We have observed the blazar Markarian 421 with the TACTIC 7-ray 
O telescope at Mt. Abu, India, from 22 November 2009 to 16 May 2010 for 265 hours. 

Detailed analysis of the data so recorded revealed presence of a TeV 7-ray signal with a 
statistical significance of 12.12(7 at i?-^ > 1 TeV. We have estimated the time averaged 
differential energy spectrum of the source in the energy range 1.0 - 16.44 TeV. The 
I spectrum fits well with the power law function of the form [dF/dE = fQE~^) with 

> /o = (1.39 ± 0.239) X 10-"cm-2s-irel/-i and T = 2.31 ± 0.14. 
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^ 1. Introduction 

o 

^ Markarian 421 falls in the blazar category of Active Galactic Nuclei (AGN) which 

are the most extreme and powerful variable sources of photons with energies ranging 

*0 from the radio to the 7-ray regimes. Blazars are believed to have their jets more 

^ aligned towards us as compared to any other class of radio loud AGN, are quite 

^ bright, have irregular amplitude variability in all accessible spectral bands and a core- 

dominated radio morphology with flat radio spectra, which join smoothly to the infra- 
red (IR), optical and ultra-violet spectra. The observed high variability and broadband 
emission from such objects make long-term multiwavelength observations of blazars 
very important for understanding their emission mechanisms and other jet related 
properties [211 162]. The Spectral Energy Distributions (SEDs) of these violent objects 
have double-peaked structure in the uF^, versus frequency plot. The first peak is usually 
referred to as the synchrotron peak, both in leptonic and hadronic models for blazar 
emission. Further, it is generally believed to be the result of incoherent synchrotron 
emission from relativistic electrons and positrons, which are conjectured to be present 
in the magnetic fields of the jet. The origin of the second peak, usually referred to as 
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the inverse-Compton peak, is less well determined. In synchrotron self-Compton (SSC) 
models [HI SHI HS], |57] it is assumed that the synchrotron photons are up-scattered to 
higher energies by the electrons while in external Compton (EC) models [HI EH l3T| IM], 
these seed photons can come from the accretion disk, broad-line region, torus, local 
infrared background, cosmic microwave background, ambient photons from the central 
accretion flow or some combination of these sources. In the framework of hadronic 
models [H |121 US], the X-ray to 7-ray emission is believed to be due to synchrotron 
radiation from protons accelerated in highly magnetized compact regions of the jet. In 
another scenario it is proposed that the proton-proton collisions, either within the jet 
itself or between the jet and ambient clouds, give rise to neutral pions which decay to 
gamma rays [Ti t [22 | Ii9]. 

The blazar Mrk 421(z=0.031) was the first extragalactic source detected at TeV 
energies in 1992 using imaging atmospheric Cherenkov telescopes [5T1 153] . The source 
has been observed extensively by various groups since then [U E] [TJ [TOl [131 [13 I23l 
[m El [IS, [521 ESI [66l [68] . These observations have shown that the TeV 7-ray emission 
from Mrk 421 is highly variable with variations of more than one order of magnitude and 
occasional flaring doubling time as short as 15 mins [SI [20]- Since its detection in the TeV 
energy range, Mrk 421 has also been the target of several multiwavelength observation 
campaigns [H [TU [T31 [13, [201 [2Z1 [SOI [SH [S3] • Several groups have also determined the 
energy spectrum of Mrk 421 at various flux levels and some of these results are cited and 
discussed in section 6. Differences in the energy spectrum of Mrk 421 and Mrk 501 have 
also been addressed to understand the 7-ray production mechanisms of these objects 
and absorption effects at the source or in the intergalactic medium due to interaction of 
7-rays with the Extragalactic Background Light (EBL) [2S]- In this paper, we present 
TACTIC results obtained from our observation campaign on this source during the 
period 2009-10. 

2. TACTIC telescope 

TACTIC (TeV Atmospheric Cherenkov Telescope with Imaging Camera) 7-ray telescope 
is located at Mt. Abu ( 24.6° N, 72.7° E, 1300m asl), Rajasthan, India. It deploys a 
F/1 type tracking light collector of ~9.5 m^ area which is made up of 34 0.6-m-diameter 
front-coated spherical glass facets which have been prealigned to produce an on-axis 
spot of ~ 0.3° diameter at the focal plane. The telescope has a 349 photomultiplier 
tubes (ETL 9083UVB) -based imaging camera with a uniform pixel resolution of ~0.3° 
and a field-of-view of ~6°x6° to record images of atmospheric Cherenkov events. Data 
used in the present work have been collected with the inner 225 pixels and the innermost 
121 pixels were used for generating the event trigger. The trigger is based on the 3NCT 
(Nearest Neighbour Non-Collinear Triplets) logic [37] . The safe anode current (< 3 
/iA) operation of the photomultipliers (PMT) has been ensured by implementing a gain 
control algorithm [15]. The data acquisition and control system of the telescope [SZ] 
have been designed around a network of PCs with the QNX (version 4.25) real-time 
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operating system. The triggered events are digitized by CAMAC based 12-bit Charge 
to Digital Converters (CDC) which have a full scale range of 600 pC. The single pixel 
threshold was set to > 14 photoelectrons(pe). The detailed description of the telescope 
related hardware and software has been given in ^9]. The telescope is sensitive to 7- 
rays above 1 TeV and can detect the Crab Nebula at 5o" significance level in 25 hours 
of observation. 

3. Mrk 421 Observations 



Table 1. Details of Mrk421 observations using TACTIC during the period 2009-10. 



Year 


Month 


Observation dates 


Total data (hrs.) 


Data Selected (hrs.) 


2009 


Nov. 


22-27 


5.6 


5.6 


2009 


Dec. 


13,15-16,18-19,21,23-27 


21.1 


16.9 


2010 


Jan. 


10,13-24 


47.8 


41.0 


2010 


Feb. 


6,10-20,22-23 


54.2 


50.2 


2010 


Mar. 


7-22 


60.5 


59.4 


2010 


Apr. 


2,4-7,9-18 


50.2 


38.5 


2010 


May 


2,4,6-13,15-16 


26.0 


18.3 


2009-2010 


Nov.-May 




265.4 


229.9 



Very High Energy(VHE) 7-ray observations of Mrk 421 presented here were made 
with the TACTIC 7-ray telescope during 22 November 2009 to 16 May 2010 for 265.4 
hours. The observations were carried out in tracking mode, where the source is tracked 
continuously without taking off- source data [5T|. This mode improves the chances of 
recording possible flaring activity from the candidate source direction. Details of these 
observations with the TACTIC 7 -ray telescope during 2009-10 are given in Table [T} We 
have used the standard data quality checks to evaluate the overall system behaviour and 
the general quality of the recorded data, including conformity of the Prompt Coincidence 
(PC) rates with the expected zenith angle trend, compatibility of the arrival times of PC 
events with the Poissonian statistics and the steady behaviour of the chance coincidence 
rate with time [65]. After applying these cuts, we have selected good quality data sets 
as per details given in Table [TJ 
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4. Data Analysis and Results 

The analysis of data recorded by the ground based atmospheric Cherenkov gamma-ray 
telescopes, involves a number of steps including filtering of night sky background light, 
accounting for the differences in the relative gains of the PMTs and finding Cherenkov 
image boundaries, image parameterization, event selection, energy reconstruction etc. 
We have characterized each Cherenkov image using a moment analysis methodology 
given in[36l [551 ESI El El E], wherein the roughly elliptical shape of the image is described 
by the LENGTH and WIDTH parameters and its location and orientation within 
the telescope field of view are given by the DISTANCE and ALPHA(q;) parameters 
respectively. In addition, the two highest amplitude signals recorded by the PMTs 
(maxl, max2) and the amount of light in the image (SIZE) are also obtained for each 
recorded image. Further, another parameter FRAC2 which is the addition of two largest 
amplitude signals recorded by PMTs divided by the image size, is also estimated [65]. 
The standard Dynamic Supercuts |3S1 US] procedure is then used to separate 7-ray like 
images from the overwhelming background of cosmic-rays. This procedure uses the 
image shape parameters LENGTH and WIDTH as a function of the image SIZE so that 
energy dependence of these parameters can be taken into account. The 7-ray selection 
criteria used in this analysis is given in Table |2] and has been obtained on the basis of 
Monte Carlo simulations carried out for the TACTIC telescope [26] 138]. 

Table 2. Dynamic Supercuts selection criteria used for analyzing the TACTIC data. 



Parameters 


Cuts Value 


LENGTH (L) 


0.11< L < (0.235 + 0.0265 xlnS)° 


WIDTH (W) 


0.065< W < (0.085 + 0.0120 xlnS)° 


DISTANCE (D) 


0.5< D < (1.27 Cos°-ssine)° (e=zenith angle) 


SIZE (S) 


S > 485 dc (6.5 digital counts = l.Ope) 


ALPHA (a) 


a < 18° 


FRAC2 (F2) 


F2 > 0.38 



VHE 7-ray signal processing is done by using the histogram of the a parameter (defined 
as the angle between the major axis of the image and the line between the image centroid 
and camera center, when the source is aligned along the optical axis of the telescope) 
after applying the set of image cuts given in Table [2j The distribution of the a parameter 
is expected to be fiat for the isotropic background of cosmic ray events, whereas for the 
7-ray signal events, the distribution is expected to show a peak at smaller a values. 
This range for the TACTIC telescope is a < 18°. The contribution of the background 
events is estimated from a reasonably fiat a region of 27° < a < 81°. The number 
of 7-ray events is then calculated by subtracting the expected number of background 
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Figure 1. Distribution of image parameter a for the TACTIC observations of Mrk 
421 during 2009-10. Vertical dashed line represents the a parameter cut used for 7-ray 
like events. Horizontal line represents the mean background level per 9° bin derived 
by using the reasonably flat a region of 27° < a < 81°. Error bars shown are for 
statistical errors only. 



events, calculated on the basis of the background region [21j, from the 7-ray domain 
events. The reason for not including the a bins 18°-27° and 81°-90° in the background 
region is to ensure that the background level is not overestimated because of a possible 
spill over of 7-ray events in the third bin and the truncation of the Cherenkov images 
recorded at the boundary of the imaging camera [21| in the last bin. The significance 
of the excess events is calculated by using the maximum likelihood ratio method of Li 
and Ma gl]. 

Next, we present data analysis results of Mrk 421 observations made using the 
TACTIC telescope during the year 2009-10. When all the data are analyzed together, 
the corresponding results obtained are shown in Figure [TJ In this figure the histogram 
of the a parameter has been shown after having applied shape and orientation related 
imaging cuts given in Table |2} It shows the presence of the VHE 7-ray signal with 
a total number of 7-ray like events of 1932.7 ± 159.4. The statistical significance of 
the detected signal stands at Vl.Vla level, thereby indicating that the source was in a 
relatively high TeV emission state during the period of observations. When these data 
were divided into seven monthly spells as depicted in Table |3] and analyzed by using 
the same data analysis procedure, the number of 7-ray like events obtained was highest 
during February 2010 observations. The details of these results have been given in Table 

El 
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Table 3. Monthly spell wise analysis of Mrk 421 data recorded during 2009-10 with 
only statistical errors. 



Spell 


Observation Period 


Time (hrs.) 


7-rays 


Significance 


1 


22-27 Nov. 2009 


5.6 


37.0 ± 20.2 


1.83 


2 


13-27 Dec. 2009 


16.9 


105.7 ± 38.4 


2.75 


3 


10-24 Jan. 2010 


41.0 


285.0 ± 62.7 


4.54 


4 


06-23 Feb. 2010 


50.2 


711.7 ± 88.6 


8.03 


5 


07-22 Mar. 2010 


59.4 


432.7 ± 88.6 


4.88 


6 


02-18 Apr. 2010 


38.5 


264.7 ± 52.9 


5.00 


7 


02-16 May 2010 


18.3 


96.0 ± 33.1 


2.90 


Total 


22 Nov. 2009-16 May 2010 


229.9 


1932.7 ± 159.4 


12.12 
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Figure 2. Mrk 421 light curves for 2009-10 TACTIC observations. 



These data have also been analyzed on nightly basis in order to look out for the 
possibility of a strong episodic TeV emission. The results so obtained are depicted in 
Figure [2k, which shows a day-to-day variations of the 7-ray rate (7-rays/hour) for Mrk 
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Figure 3. Distribution of image parameter a for the TACTIC observations of Mrk 
421 during February 2010. Vertical dashed line represents the a parameter cut used 
for 7-ray like events. Horizontal line represents the mean background level per 9° bin 
derived by using the reasonably flat a region of 27° < a < 81°. Error bars shown are 
for statistical errors only. 

421 during 2009-10 observations. This light curve is characterised with a reduced 
value of 83.26/69 with respect to the mean level of 7.4±0.65 photons per hour, with 
corresponding probability of 0.12 which is consistent with the no variability hypothesis. 

Next we compare the TeV light curve obtained which is shown in Figure [2^ 
with those of the source with the Swift /'Quisi Alert Telescope (BAT) |5S|( 15-50 
KeV), Rossi X-ray Timing Explorer (RXTE) /All-Sky Monitor (ASM) [12] ( 2-10 KeV) 
and Fermi /Large Area Telescope (LAT) [28] (1-300 GeV) detectors. The Swift /BAT 
contemporary light curve of the source obtained from its archived data [5H] is shown in 
Figure [2}3. It is characterised with a reduced value of 1594/168 (probability is close 
to zero which is consistent with the variability hypothesis) with respect to the mean 
level of 0.00349 ± 7.4e-05 counts cm-^g-^ Next the RXTE/ASM contemporary light 
curve shown in Figure ^ which has been plotted by using the daily average count rates 
of the ASM from its archived data [12]. The ASM light curve is characterised with a 
reduced value of 2825/174 ( probability obtained is very low which is consistent with 
the variability hypothesis) with respect to the mean level of 2.15 ± 0.025 counts/sec. 

The Fermi /LAT contemporary light curve shown in Figure |2]i has been plotted by 
using the daily counts cm~^s~^ of the LAT from its archived data [2H]- The LAT light 
curve is characterised with a reduced value of 291/156 ( probability obtained is very 
low which is consistent with the variability hypothesis) with respect to the mean level 
of 6.17e-08 ± 1.25e-9 counts cm~'^sec~^ . 
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Figure 4. Mrk 421 light curves for February 2010 observations. 



As already mentioned above and is also clear from Table [3] that we have detected 
a statistically significant TeV 7-ray signal at 8.03(7 level during February TACTIC 
observations of the source. It was observed for 12 nights and a total of 50.2 hours 
of on-source observations were possible during this month. Total number of TeV 7-rays 
detected during this month are 711. 7± 88.6 which is clear from the corresponding alpha 
plot shown in Figure [3] Further, we find that the source was at higher TeV emission 
state during three nights of 15- 17 February, 2010(MJD 55242.73913-55244.98546). 
Details of these results are given in Table |4} On 16th February, 2010 we detected 
177.7 ± 30.8 7-ray photons of TeV energies at more than 5cr statistical significance 
level. During other two nights of 15th and 17th February, 2010 again we detected 
statistically significant number of TeV 7-rays with greater than 4(T level. Using these 
results we believe that the source had possibly switched to an intense fiaring state 
during 15- 17 February, 2010 observations. Further, in order to compare the February 
2010 TeV results with those at lower energies, we have depicted the February source light 
curves of TACTIC, Swift /BkT, RXTE/ASM and Fermz/LAT in Figures |4|i,b,c and d 
respectively. As is clear from these figures the daily variations in the number of counts 
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Table 4. Details of Mrk 421 data analysis results of TACTIC observations made 
during 15-17 February 2010. 



Obs. Date 
(Feb. 2010) 


Start MJD 


End MJD 


7-rays 


Sig.(a) 


Obs. Time(hrs) 


15 


55242.73913 


55242.99233 


93.3 ± 27.5 


3.39 


4.9 


16 


55243.77758 


55243.98686 


177.7 ± 30.8 


5.76 


5.0 


17 


55244.75948 


55244.98546 


143.0 ± 30.2 


4.73 


5.2 


15, 16, 17 


55242.73913 


55244.98546 


414.0 ± 51.2 


8.09 


15.1 



detected by these four experiments are such that when we fit zero degree polynomials to 
these hght curves the reduced values obtained are 33.46/11, 353.8/15, 484.1/14 and 
22.38/14 (xVdof ) for TACTIC, Swift /BAT, RXTE/ASM and Fermi /hKT respectively. 
Fermz/LAT's light curve reduced value of 22.38/14 yields a corresponding probability 
of 0.07 thereby indicating that the hght curve is consistent with a fiat distribution. 
However, the first three light curves yield corresponding probabilities much less than 
0.05 thereby indicating in favour of variable signal detection in the TACTIC, Swift /BAT 
and RXTE/ASM energy ranges. These light curves also indicate a possible correlation 
between TeV and X-ray emissions from this source direction particularly during 15- 17 
February, 2010 observations. This lends possible support to the SSC model in which a 
unique electron population produces the X-rays by synchrotron radiation and the 7-ray 
component by inverse Compton scattering. 

Further when hourly data runs of these three days are investigated in terms of the 
detected number of 7-ray like events, we find that on the night of 17th February, 2010, 
during the hourly run from MJD 55244.812502 to 55244.854161, we have detected 57 ± 
14.2 7-ray photons of TeV energies at 4a statistical significance as is clear from Figure 
|5] It may be noted here that TACTIC detects on an average 9 7-ray like events in one 
hour of observations from the Crab nebula direction. Therefore, we conclude that the 
source has possibly exhibited yet another short duration high TeV 7-ray emission state 
during this period. 

5. Differential energy spectrum 

In this section we present the differential time-average energy spectrum of the source 
estimated using its TACTIC observations made during the period 22 Nov. 2009- 
16 May 2010. The energy reconstruction procedure employed to unfold it uses an 
Artificial Neural Network (ANN) with 3:30:1 configuration and resilient backpropagation 
algorithm. The energy of each 7-ray like event is estimated on the basis of its image 
SIZE, DISTANCE and zenith angle. This method yields an energy resolution of ~ 26%, 
the details of which have been given in lUB] . 

The observed Mrk 421 differential energy spectrum so obtained has been tabulated 
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Figure 5. TACTIC observed hourly light curves of Mrk 421 on three nights of 15,16 
and 17 February 2010 (MJD 55242.73913 to 55244.98546). 



Table 5. Differential energy spectrum data for Mrk 421 in 2009-10 with the TACTIC 
telescope. Only statistical errors are given below. 



Energy 

(TeV) 


Diff. flux 

photons cm^^ s^^TeV^^ 


Statistical error in flux 

photons cm-2 s-i TeV-i 


1.00 


1.75 X 10-11 


6.78 X 10-12 


1.50 


5.45 X 10-12 


1.11 xlO-12 


2.22 


2.15 xlO-12 


3.15 xlO-13 


3.32 


7.46 X 10-13 


1.54 xlO-13 


4.95 


3.89 X 10-13 


7.98 xlO-i^ 


7.38 


1.85 X 10-13 


4.31 xlO-i^ 


11.00 


5.11 xlO-i^ 


2.59 xlO-i^ 


16.44 


1.01 xlO-i^ 


1.24 xlO-i^ 



in Table |5] and also shown in the Figure [6] which extends up to 16.44 TeV. The errors in 
the flux constant and spectral index are standard errors. A power law function of the 
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Figure 6. TACTIC observed differential energy spectrum of Mrk 421 derived using 
the 2009-10 observations( filled circle). De-absorbed source spectra estimated by using 
various EBL models are also shown for comparison. 



form dF/dE = f^E ^ when fitted to the data in the energy range 1- 16.44 TeV yields 

^ = (1.39 ± 0.239) X lO-iicm-2s-irer-i(^-)-2-3i±o.i4 
dE l.OTeV 

for the observed differential energy spectrum of the source with x^/dof = 3.37/6 

(probability 0.76). 



Table 6. Estimated intrinsic differential energy spectrum fitted parameters for three 
models. 



r 


Ar 


fo 


Afo 


Model 


2.31 


0.14 


1.39 X 10-^^ 


2.39 X 10-12 


Presented observations 


2.03 


0.16 


1.86 X 10^^^ 


3.60 X 10-12 


Primack (50j 


1.85 


0.16 


1.73 X 10-11 


3.16 X 10-12 


Stecker Baseline [59] 


1.73 


0.17 


2.00 X 10-11 


4.11 X 10-12 


Stecker Fast evolution [59] 


1.94 


0.18 


1.60 xlO-ii 


3.43 X 10-12 


Franceschini [29] 


1.97 


0.14 


1.88 X 10-11 


3.20 X 10-12 


Gilmore [32j 



We have also estimated the source intrinsic 7-ray spectra by using the SED of the 
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Figure 7. Comparison of the observed spectral energy distribution of Mrk 421 in the 
VHE range obtained by various groups. TACTIC observed differential energy spectrum 
during 2009-10 observations is in the energy range 1.0 - 16.44 TeV. 



EBL provided by five models given in Table |6} The resulting intrinsic spectra for five 
models are shown in Figure [6j A power law fit to the estimated intrinsic data of the 
type d^/dE = /o^"^ in the energy range 1- 16.44 TeV yields values of T and /o listed 
in Table [6j As is clear from Figure |6} the Primack model yields minimum attenuation 
of TeV photons, in contrast, Stecker's fast evolution model yields maximum attenuation 
as compared to other models used in estimating the intrinsic source spectra. 

As mentioned earlier, Mrk 421 is a nearby source with z=0.031, so we do not 
expect significant absorption of intrinsic source spectrum due to EBL in the TACTIC 
energy range. Another important characteristic feature associated with this source 
is its relatively higher flaring frequency when compared to other objects of the same 
class, thereby making it an excellent cosmic laboratory to study VHE 7-ray emission 
mechanisms and physics of relativistic jets. 

6. Discussion and conclusions 



For the reasons mentioned above, we are monitoring Mrk 421 on the long term basis 
with TACTIC 7-ray telescope. In this paper, we have presented the results obtained 
using TACTIC observations made during 2009-10. Detailed analysis of this data shows 



p. Chandra et al. 



13 



evidence for the presence of a statistically significant VHE 7-ray signal at 12.12(j level. 
The total number of 7-ray like events detected are 1932.7 ± 159.4. These observations 
suggest that the source had gone into a variable active state in the VHE region during 
the period of TACTIC 2009-10 observations. During the month of February 2010, TeV 
signal strength in terms of the number of TeV 7-ray like events which is 711. 7± 88.6 
highest as compared to results obtained for other months as depicted in Table |3] It is 
interesting here to mention that the HESS [S], VERITAS |17] and ARGO-YBJ [HIES] 
have also detected the source in a high state particularly during February 2010. 

We have estimated an observed differential energy spectrum of the detected 7-ray 
like events which is shown in Figure [6j The spectrum fits well with the power law 
function of the form {dF/dE = /o^"^) with /o = (1.39 ± 0.239) x 10'^^ cm-^s-^TeV-^ 
and r = 2.31 it 0.14. In comparison with our earlier TACTIC detected high states of 
this source reported in |66l EH], present observed spectrum is harder and extends up 
to 16 TeV primary 7-ray photon energy. We have also used five templates of EBL to 
estimate r's for z = 0.031 which have been used to estimate corresponding intrinsic 
source spectra. The results of this study have been tabulated in Table [6] 

A comparison of the observed spectral energy distribution of this source in the VHE 
range obtained by the HESS, MAGIC, VERITAS, Whipple, HEGRA, CAT, ARGO-YBJ 
and TACTIC [HIiniEHlllE] [SlE21ESlE31lll]is shown in Figure (T) As is clear from 
this figure our 2009-10 results are in close agreement with various reported results on 
this source. 

TACTIC, Swzft{BAT) and RXTE(ASM) light curves of 15- 17 February, 2010 
observations (55242.73913 - 55244.98546 ) do indicate a possible correlation of TeV 
and X-ray emissions of the source as is shown in Figures |4^,b and c. This indicates 
that the mentioned short TeV 7-ray flare was not an orphan type as it is also seen at 
lower energies and possibly supports the SSC model in which a relativistic electrons 
population produces the X-rays by synchrotron radiation and the 7-ray component by 
inverse Compton scattering. 
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